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f*^) ■ Abstract. We investigate the properties of highly compressible turbulence and its 

C*~) | ability to produce self-gravitating structures. The compressibility is parameterized 

by an effective polytropic exponent j c . In the limit of small 7 e , the density jump 
at shocks is shown to be of the order of e M , and the production of vorticity by the 
nonlinear terms appears to be negligible. In the presence of self-gravity, we suggest 
| that turbulence can produce bound structures for j c < 2(1 — 1/"-), where n is the 

£Sl ■ typical dimensionality of the turbulent compressions. We show, by means of numerical 

simulations, that, for sufficiently small j e , small-scale turbulent density fluctuations 
eventually collapse even though the medium is globally stable. This result is preserved 
in the presence of a magnetic field for supercritical mass-to-flux ratios. 
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& : INTRODUCTION 



In this paper we present a brief discussion on the dynamical properties of highly 
compressible turbulence, with and without self-gravity. A full-length discussion can 
be found in Vazquez-Semadeni et al. (1996, hereafter Paper I). The compressibil- 
ity of the medium is parameterized by an effective polytropic exponent 7 e arising 
from balance between heating and cooling processes (Elmegreen 1991; Vazquez- 
Semadeni et al. 1995; Passot et al. 1995), such that the pressure P and the density 
p are related by P oc p 7c . We consider the cases where 7 e is either constant, or has 
a piecewise density dependence (a "piecewise polytropic model" or ppm). We also 
consider fully thermodynamic cases. Most of the results are based on numerical 
calculations which solve the full MHD equations in two or three dimensions, at 
resolutions of 128 2 or 64 3 , respectively (Paper I). 
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NO SELF-GRAVITY 

Density Structures 

It can be easily shown that the density jump X = P2/P1 across a shock in a 
barotropic gas of index 7e satisfies: 

X 1+ >-(l+ 7c M 2 )X + 7c M 2 = 0, (1) 

where M is the upstream Mach number. 

In the isothermal case where j e — 1, X — M 2 , while for < 7e <C 1, X ~ e M . 
This implies that for small 7e : 

i) The density jump is much larger than in the isothermal case. 

ii) Less than supersonic motions (with respect to the isothermal sound speed) are 
sufficient for producing large density fluctuations. 



Vorticity evolution 

From the momentum conservation equation (Paper I), we can derive the evolution 
equation for the potential vorticity u p = u/p = V x u/p : 

^ + u • Vcj p = uo p • Vu + V x F s + VP x Vp/p 3 (2) 

where F s is the solenoidal (or rotational) part of the turbulent energy sources. 

While compressive motions can easily be generated from vortical motions, only 
the "vortex stretching" term u p ■ Vu and the "baroclinic" term (nonzero only in the 
fully thermodynamic case) are available as nonlinear sources of u p . How effective 
are they? 

Fig. 1 (left) shows the evolution of the ratio e s /ek of solenoidal to total kinetic en- 
ergy per unit mass for four 3D runs. Runs 60 and 67 are ppm runs with 7em i n = 0.25 
and 7 e m i n = 0.75, both with fully solenoidal initial velocity modes but fully com- 
pressible forcing (F s = 0). Run 68 is similar to run 60, but fully thermodynamic. 
Run 69 is similar to run 60, but with a Coriolis term added. In all cases except run 
69, e s /ek is seen to decay at roughly the same rate in time. Thus : 

i) We have found negligible nonlinear transfer from compressible to rotational ki- 
netic energies. Similar results in the weakly compressible case have been found by 
Kida & Orszag (1991a,b). 

ii) Additional sources of vorticity, such as the Coriolis force (acting on large scales), 
are necessary for the maintenance of significant amounts of vorticity. 

Additionally, the presence of an initially uniform magnetic field also allows for the 
production of rotational energy from the compressive motions, with an efficiency 
proportional to the magnetic field strength (fig. 1 (right) ; in runs 86, 84, 87, 88, 
the non-dimensional field intensity is respectively 0.05, 0.3, 1 and 3). This process 
is more important at small scales. 
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FIGURE 1. Evolution of the ratio of rotational to total specific kinetic energy for various runs. 
Left: non-magnetic runs. Right: magnetic runs. 

FLOWS WITH SELF-GRAVITY 

Non- magnetic case 

It can readily be shown that for a barotropic medium with polytropic exponent 
7e, the effective Jeans length is : 



^eff 



1/2 




(3) 



where q is the isothermal sound speed, such that P = cfp, and Lj is the Jeans 
length based on q. 

The critical density required for destabilizing a length scale L is thus pj oc 
£2/(7 -2)^ j n orc [ er t account for turbulent compressions acting on n directions, we 
consider a volume V = L n Ll~ n , where L is the side of the volume which varies upon 
compression, and Lq is the side that remains unaltered. The critical mass to desta- 

ra+2 „_ 

bilize this volume is thus given by Mj oc L^- 2 Lq n . If Mj is a decreasing function 
of L, turbulent compression can produce gravitationally unstable structures, i.e, if 
(see also McKee et al. 1993) 



7c < 7crit = 2(1 - l/n). 



(4) 



This result is illustrated in fig. 2 (left), which shows the evolution of the global 
density maximum for three single-7 e simulations, all with an isothermal Jeans 
length Lj = 1.1 times the box size (i.e., gravitationally stable) and purely com- 
pressible random forcing, but with 7 e = 0.9, 0.3 and 0.1 respectively. While the 
simulation with 7 C = 0.9 never develops a gravitational collapse, the other two do, 
earlier in the case with 7 e = 0.1. 
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FIGURE 2. Left: Evolution of the global density maximum for three runs with 7 C = 0.9 (solid 
line), 7o = 0.3 (dotted line) and 7 e = 0.1 (dashed line). Right: Final density fields of two magnetic 
simulations, one with j e = 0.9 (center) and one with j e = 0.3 (right). 

Magnetic case 

The above result is preserved in the presence of a magnetic field. Fig. 2 also 
shows contour plots of the final density field in simulations with an initially uni- 
form magnetic field along the x-direction, a Jeans length Lj = 0.9 (i.e., gravi- 
tationally unstable with respect to a uniform density configuration), and purely 
compressible random forcing. One simulation has j e = 0.9 (center) and the other 
has 7 e = 0.3 (right). The latter is seen to have undergone collapse, while the former 
only contracts to a pancake-type structure supported by thermal pressure against 
final collapse. However, for large enough magnetic strengths, we have found that 
turbulence is again not capable of inducing gravitational collapse, suggesting that 
a subcritical magnetic regime (Mouschovias & Spitzer 1976) cannot be forced to 
collapse by large-scale turbulence (although the presence of small-scale turbulent 
modes may affect this result). 
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